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Experimental observations on synaptic plasticity at individual glutamatergic synapses from the 
CAS Shaffer collateral pathway onto CAl pyramidal cells in the hippocampus suggest that the 
transitions in synaptic strength occur among discrete levels at individual synapses ( |^ and S. 
S.-H. Wang, unpublished data used with the authors' permission). This happens for both long term 
potentiation (LTP) and long term depression (LTD) induction protocols. O'Connor, Wittenberg, 
and Wang have argued that three states would account for their observations on individual synapses 
in the CA3-CA1 pathway. We develop a quantitative model of this three state system with transi- 
tions among the states determined by a competition between kinases and phosphatases shown by 
O'Connor et al. to be determinant of LTP and LTD, respectively. Specific predictions for various 
plasticity protocols are given by coupling this description of discrete synaptic AMPA conductance 
changes to a model of postsynaptic membrane potential and associated intracellular calcium fluxes 
to yield the transition rates among the states. We then present various LTP and LTD induction 
protocols to the model system and report the resulting whole cell changes in AMPA conductance. 
We also examine the effect of our discrete state synaptic plasticity model on the synchronization of 
realistic oscillating neurons. We show that one-to-one synchronization is enhanced by the plasticity 
we discuss here and the presynaptic and postsynaptic oscillations are in phase. Synaptic strength 
saturates naturally in this model and does not require artificial upper or lower cutoffs, in contrast 
to earlier models of plasticity. 
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I. INTRODUCTION 



Experiments on synaptic plasticity at individual 
synapses in CAS to CAl hippocampal pathways reveal 
an "all-or-none" change in their synaptic strength ( [2ll| 
and S. S.-H. Wang, unpublished data used with the au- 
thors' permission). Indications of this were seen a decade 
ago (C. F. Stevens and Y. Wang, personal communica- 
tion). The recent measurements have given substantial 
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standing to the notion that single synapses may operate 
as a discrete state system in their plastic changes asso- 
ciated with long-term potentiation (LTP) and long-term 
depression (LTD). 

In this paper we first explore a general formulation of 
the possibility that individual synapses can express a fi- 
nite number of discrete levels of conductance rather than 
the continuous or g raded or analog picture often formu- 
lated P, 0, 0, l25| . [2l| have commented on the posi- 
tive consequences for reliability of neural memory from 
discrete state synapses. In addition there is computa- 
tional evidence that at an individual dendritic spine the 
number of free Ca^+ions at the peak of a Ca^+transient 
evoked by synaptic stimulation may be small jl2j . A 
synapse which must respond to this very small signal 
might well select a strategy of an "all-or-none" response 
in adjusting its strength as a means of achieving a mea- 
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sure of reliability. 

After a brief discussion of a general formulation of an 
L level synapse, we focus our attention on L = 3, which 
is suggested by the recent data of O'Connor et al. (D. 
H. O'Connor, G. M. Wittenberg, and S. S.-H. Wang, un- 
published data used with the authors' permission). With 
three levels we explore the transitions among the levels 
using observations of O'Connor et al. With some gen- 
eral arguments based on the measurements, we are able 
to establish values for the normalized conductances of 
the individual L = 3 levels, and we suggest an exper- 
iment which will determine an interesting ratio among 
the transition rates. 

The transition rates themselves depend on a model 
for how complex chemical pathways lead to the change 
in AMPA conductance and the number of AMPA recep- 
tors at any synapse H S [ll Ellll 113 ■ We adopt a 
version of our earlier model 0] of these processes to pro- 
vide a basis on which to make quantitative predictions 
of the outcome of various LTP /LTD induction protocols 
on the changes in AMPA conductance of a postsynap- 
tic cell. We explore spike time dependent protocols as 
well as presentation of spike bursts of various frequencies 
to the postsynaptic cell. These compare well to experi- 
ment, in particular to results presented by |23, and pre- 
dictions are made for various spike timing experiments. 
We also explore the following setup: a periodically fir- 
ing conductance based presynaptic neuron provides ex- 
citatory synaptic input to a periodically firing conduc- 
tance based postsynaptic neuron. In certain ranges of 
frequency and conductance strength for the excitatory 
connection these neurons synchronize [22|. We explore 
the effect on this synchronization of our discrete synaptic 
strength model and show that the regime of one-to-one 
synchronization is significantly enlarged |32| for some 
conductance strengths and that the synchronization is 
an in-phase firing of the two neurons. 

It is an important feature of models built on a fi- 
nite number of discrete synaptic levels that the synaptic 
strength (AMPA conductance) is always bounded above 
and below. This is in contrast to models developed over 
the years, including our own pol |- which do not have 
this property. 



II. METHODS 



A. Discrete State Synapses; Transition Rate 
Models 



The data of O'Connor et al. suggest that three dis- 
crete states of AMPA conductance are found at individ- 
ual synapses. There are suggestions of other numbers of 



discrete levels ( |2l|; C. F. Stevens and Y. Wang, per- 
sonal communication). If the total number of levels is 
L and is indexed by ^ = 0, 1, 2, L — 1, and there are 
Ns synapses indexed by n = l,2,...,iVs, then we rep- 
resent the occupation of synapse n in state 1 at time t 
by Np (t) . These occupation numbers are either zero or 
unity. They can change in time due to LTP /LTD induc- 
tion protocols or other biological processes. The average 
occupation number in state I is given by 



1 Ns 



it). 



(1) 



These pi(t) will constitute the main dynamical vari- 
ables of our model. They are taken to satisfy linear rate 
equations of the form 



dpijt) 
dt 



fc=0 ^ 



The transition rates {t), Wi^i^k) (t) are selected to 

assure 



L-l 



(3) 



fe=0 



As long as the number of active synapses Ng is un- 
changed, this last equation follows from the definition 
oi Pk{t)- In the limit of Ns large, we assume that the 
Pi{t) remain finite; this is a standard assumption about 
such a limit for the description of a large number of ob- 
jects each having a discrete set of states. 

One can collect the average occupation numbers into 
an L-dimensional vector P(t) = (jio{t),pi{t), ...pL-iit)) 
which satisfies 



dP{t) 
dt 



MW-P(i), 



(4) 



and the conservation rule means that M(t) always has at 
least one zero eigenvalue 14, 28j the dynamics of P(i) 
takes place in the L — 1 dimensional space orthogonal to 
the constant L-dimensional vector C = (1,1,1,...,!). 

The effect of having this constraint may be seen in the 
decomposition of P(i) = C-|-P_L(i) where P_L(t) • C = 0. 
The dynamics of motion for Pj_ (t) is 



dP±{t) 
dt 



M{t)-F^{t)+M{t)-C, 



(5) 



which is driven motion of the vector P±{t) spanning the 
L — 1 dimensional space orthogonal to C. P±(i) is de- 
fined up to a rotation about C. This dynamical descrip- 
tion is similar to that of driven precession of a spin in a 
time dependent magnetic field. 
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Each discrete level I — 0, — 1 lias an AMPA 

conductance g; normalized to baseline. The normalized, 
diniensionless AMPA conductance of the neuron with 
Ns synapses is 

Ns L-1 

GAMPAit) = Y.9iNrit) 

n=l 1=0 
L-1 

= NsY,9iPi{t). (6) 



This means that the quantity 



Ga 



MPA 



it) 



L-1 



Ns 



(7) 



1=0 



the normalized, dimensionless AMPA conductance per 
synapse, is independent of the number of synapses when 
Ng is large, depending only on the average occupa- 
tion number of the synaptic levels and the conductance 
associated with each level. By definition, before any 
LTP /LTD induction protocols are presented to the post- 
synaptic neuron, this quantity is equal to one: 



GAMPAjt = 0) 
n's 



1. 



(8) 



i(t) 



Ns 



1 



In our work below, we report the quantity 
as the output from our simulation of various induction 
protocols. 

To fully specify the model of synaptic plasticity as- 
sociated with the presence of discrete levels, we must 
identify the conductances gi of each level, and, through 
some form of dynamics of the postsynaptic neuron, de- 
termine the transition rates. In the next section we do 
this for the suggested three state model of O'Connor 
et al. However, if observations indicate that there are 
L ^ 3 levels operating at some synapses, then the gen- 
eral formulation presented here will cover that situation 
as well. 



B. Three State Model 

If there are three states I = 0, 1,2 then we need to 
identify three discrete level conductances go,gi,g2 and 
the transition rates among the levels. O'Connor et al. 
call their three levels "low" (state here), high (state 1 
here), and "high locked- in" (state 2 here). They suggest 
that transitions associated with LTD protocols connect 
state 1 to state 0, and transitions associated with higher 
frequency protocols, typically leading to LTP, connect 
state to state 1 and state 1 to state 2. They also 



note that when an LTD protocol is applied following a 
saturating LTP protocol to a population of synapses, the 
synapses cannot be depressed as fully as when the LTD 
protocol is applied to a naive population of synapses. 
The amount of depotentiation possible decreases over 
the 10 minutes following LTP induction. This led them 
to suggest the presence of a "high locked-in" state, called 
H* by them; we call this state 2. They do not require 
transitions between state 2 and state to account for 
their data, and we assume none as well. 

Using these observations we associate an "LTD tran- 
sition rate" g{t) with the transition between state 1 
and state 0. Similarly we associate an "LTP transition 
rate" f(t) with transitions between state and state 1. 
Loosely speaking we think of g{t) as an aggregated ac- 
tion of phosphatases leading to the dephosphorylation 
or removal of synaptic AMPA receptors and f{t) as the 
aggregated action of kinases operating in the opposite 
fashion 0, |33| . In the next section we will specify how 
one evaluates these transition rates from a dynamical 
model of the postsynaptic neuron, but for the moment 
we note that / and g will depend on the elevation of 
intracellular postsynaptic calcium concentrations above 
the equilibrium level Co ~ 100 nM. Denoting the time 
course of postsynaptic intracellular calcium concentra- 
tion as Ca(t), we define 



ACa{t) 



Ca{t) - Co 
C'o ■ 



(9) 



and the transition rates f{t),g{t) are determined by 
ACa{t) in a manner specified in the next section. 

The transition between state 1 and the "high locked- 
in" state called 2 is taken to be proportional to the 
0^1 transition rate f{t). If one had more detailed 
information on the biophysical kinase and phosphatase 
pathways, one could replace this simple assumption by 
a more complex quantity. We take this transition rate 
as bf{t) with b a constant to be determined. 

The picture outlined by O'Connor et al. does not sug- 
gest a transition from state 2 to state 1, but we find it 
is necessary. For the moment we call this transition rate 
h{t), and we will argue that it is proportional to f{t). 
h{t) cannot be zero, if the transition rate framework is 
to be consistent with observations. 

This discussion leads us to the transition rate (or 
"master" ) equations associated with the scheme depicted 
in Figure 1: 



dpojt) 

dt 
dpijt) 

dt 
dp2{t) 

dt 



^f{t)po{t)+g{t)pi{t) 

f{t)po{t) + h{t)p2{t) - g{t)pi{t) - bf{t)pi{t) 

bf{t)p,{t)-h{t)p2{t). (10) 
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By construction 

dipoit) 



- Pl{t) + P2{t)) 



dt 



0. 



(11) 



Under prolonged stimulation the postsynaptic intra- 
cellular calcium levels reach approximately constant val- 
ues, and we can ask what is the behavior of P(i) = 
[po{t),pi(t),p2{t)] under such circumstances. This 
means the functions f{t) and g{t) are thought of now 
as constant in time with magnitude determined by the 
saturated level of Ca^+. We associate this value of P, 
after long constant Ca^+elevation, with the fixed point 
of the equations (|10|l for P{t). If the induction protocol 
lasts a time T/, then the state long after T/ starting with 
an initial state P(0) will be the fixed point 



fixed point 



{9hjh,bp) 
Hf + 9)+ bp 



(12) 



In our models, low values of saturated intracellular cal- 
cium elevation ACa are connected with LTD and higher 
values, with a competition between LTD and LTP . 
The specific form of the connection between ACa and / 
and g will be given shortly, but their general dependence 
is shown in Figure 2 Q • In an LTD protocol g ^ but 
/ « 0. In an LTP protocol both / and g may be nonzero. 

If we take h = ag, then a saturating LTD protocol, 
with g 7^ 0, / w 0, will deplete both states 1 and 2, 
leading to a final state P = (1,0,0). 

If, however, we apply a saturating LTP protocol where 
neither / nor g is zero, thus arriving at (|12|l . and then 
apply a saturating LTD protocol, the choice h = ag will 
lead us back to the state P = (1, 0, 0), which is not what 
is observed. Indeed, O'Connor et al. note that the state 
reached by a saturating LTP protocol depotentiates to 
a state intermediate between all synapses in state and 
the fully saturated state; namely, our fixed point ((T^ . 

If we choose h = af, then this saturating LTD protocol 
following the saturating LTP protocol leads us to 



(a(/-fg),0,b/) 
«(/ + <?) + 6/ ' 



(13) 



namely we depopulate state 1 due to the action of g. 
This is the kind of depotentiated, but not baseline, state 
seen by O'Connor et al. 

We conclude that the choice h = af is consistent with 
the observations, and we cannot have a = 0. If a = 0, 
the high locked state would be totally populated by a 
strong LTP protocol and the synapse would not leave 
that state. Indeed, O'Connor et al. indicate that after 
such a strong LTP protocol (two rounds of theta-burst 
stimulation), most but not all synapses, about 80%, are 
in state 2. 



This completes the general formulation of the three 
level transition rate model. We now turn to the de- 
termination of the AMPA conductances gi in each level 
/ = 0, 1, L — 1 from the data presented by O'Connor et 
al. Then we discuss a conductance based neural model 
for the postsynaptic cell which permits us to translate 
electrophysiological activity into transition rates useful 
in the equations of P(t). 



C. Numerical Method 

All the simulations for the model presented in this 
work were written in C and used a 4th order Runge- 
Kutta algorithm with a fixed time step of 0.01 ms. They 
were run under Linux on a computer with an Athlon 
2400 MHz processor. 



III. RESULTS 

We first establish, using the measurements of 
O'Connor et al., the values of the normalized, dimen- 
sionless AMPA conductances of the three levels at an 
individual synapse. Our arguments show that they are 
determined independently of the specific model for the 
transition rates. Then we develop a model for the transi- 
tion rates which will allow us to make predictions about 
the response of the cells to various LTP/LTD protocols. 



Determination of the Discrete Level 
Conductances 



At t = the observed average occupation of levels is 
observed to be P(0) — (|, j, 0). This means the normal- 
ized AMPA conductance is 



G 



AMPA 



(0) 



Ns 



3go 

4 



3i 
4 ■ 



(14) 



If a strong, saturating LTD protocol is applied to 
this state, we reach P — (1,0,0), where the normal- 
ized AMPA conductance is go. It has been observed 
(O'Connor et al.) that after the induction ^^^^ = 

0.65 ± 0.03. We take this to be ^ 



Ns 



go which 

implies gi = 2. 

Next apply a phosphatase blocker (okadaic acid was 
used by O'Connor et al.) so g = and, as they did, 
present a saturating LTP signal to arrive at the state 

^ ~ ^(a+b)/^ ~ ''{a+b) ' '^'^i*^'^ independent of the tran- 
sition rate /. This is precisely the fixed point noted 
above with g — and h = af. After this protocol. 
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the normalized AMPA conductance is approximately 2, 
leading to 



agi + bg2 



= 2, 



(15) 



and thus g2 = 2. 



normalized individual level conductances {go — g 



Our model corresponds to the set of 

' 51- 

2,52 = 2). 

The constants a and b are not determined by the ob- 
servations so far. Presumably they can be determined 
by applying various induction protocols once we have a 
model for the transition rates. The actual time series of 



^^^^^lMt) + 2Mt) + 2p2{t) (16) 

will depend on a and b. O'Connor et al. provide some 
evidence that the ratio |- could be about j, but it is 
not clear saturating protocols were used in these exper- 
iments. 

This ratio ^ can be determined by another experi- 
ment not yet conducted. Start with the naive synapse 
= apply okadaic acid, so g = 0, (as 

in O'Connor et al.) which blocks phosphatases, and 
present a saturating LTP protocol. This leads to the 
state P = . Now wash out the okadaic acid and ap- 

ply the kinase blocker k252a, setting f — 0, and present 
a saturating LTD protocol. This leads one to the state 
. The normalized AMPA conductance in this state 

is 



Gampa ago + bgi 



Ns 



3 b ^ ^ 
1 + * ■ 



(17) 



after this protocol sequence would 
|. If I were 0.25, as suggested above, 

1.73. 



Measuring i^^p^ 
give us a value for 

then this normalized conductance would be || 
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B. Transition Rates in a Model for Voltage 
and Calcium Dynamics 

Evaluation of the transition rates / and g requires a 
specific model describing how the postsynaptic cell re- 
sponds to various induction protocols presented either 
presynaptically or as paired presynaptic and postsynap- 
tic actions. It also requires a model for the dynam- 
ics of Ca^+in the postsynaptic cell. We proceed using 
the idea that changes in AMPA conductance are in- 
duced by the time course of elevation of intracellular 
Ca2+ H El m m I23. The details of the pathways 
which follow elevation of Ca^+are not specified in the 
phenomenological approach we use. 



In this regard we have explored both one and two com- 
partment models of the voltage and intracellular cal- 
cium dynamics of a cell with AMPA receptors whose 
strength is changed as a result of biochemical pathways 
activated by the induction protocols. The two compart- 
ment model, which we utilize here, separates the cel- 
lular dynamics into a somatic compartment where ac- 
tion potentials are generated by the familiar sodium and 
potassium currents and a dendritic spine compartment 
where AMPA and NMDA receptors are located and in- 
tracellular calcium dynamics occurs. The details of this 
model are located in the Appendix to this paper. As one 
improves this model or replaces it with further experi- 
mental insights into the processes involved in LTP /LTD 
induction, one can use those improvements to provide 
evaluations for the transitions rates needed in the dis- 
crete state synapse model. 

The output from the biophysical model of the neuron 
which we require in this section is focused on the time 
course of elevation of intracellular calcium concentration 
relative to the equilibrium concentration Co « 100 nM. 
We call this time course Ca(t) 
way in which 



[Ca'^~^]i{t) and seek the 



ACa{t) 



Cajt) - Co 

C'n 



(18) 



influences the transition rates f{t),g{t). 

Our model involves two auxiliary variables, P(t) and 
D{t), which satisfy first order kinetics driven by Hill 
functions dependent on ACa{t). These variables satisfy 



dP{t) 

dt 
dD{t) 

dt 



Fp{ACa{t))il ^ Pit)) 



FDiACa{t))il - Dit)) 



m 

Tp 

m 



TD 



with driving terms 
Fpix) 



Fd{x) 



ax 



M 



cM 



(19) 



(20) 



We used the constants Tp — 10 ms, td = 30 ms, a = 
1.5, L = 8,M = 4.75, = 5.5, and = 10.5 in our 
calculations for this work. These equations are discussed 
in our earlier work 

These kinetic quantities are driven by elevation in 
Ca^+, ACa{t) > 0, from their resting value of zero. They 
are taken to be related to the transition rates as 



/(i,ACa(0) = P{t)D{t)^ 
g{t,ACa{t)) = P{t)'^D{t), 



(21) 



and 77 = 4 as used in our earlier work. The quantities 
f{t) and g{t) have dimensions of frequency. Our argu- 
ments do not establish their magnitude but only provide 
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a connection to their dependence on elevation of intracel- 
lular Ca^+levels. Multiplying the relations here between 
f{t) and g{t) and P{t) and D{t) by a constant rescales 
the time while not affecting the final states which lead to 
specific statements of AMPA conductance changes after 
an induction protocol. 

The model for voltage dynamics and Ca^+dynamics 
is now established. To proceed we specify an electro- 
physiological protocol. For example we present a burst 
of spikes to the presynaptic terminal with an average 
interspike interval (ISI) of our choice. Our presynap- 
tic terminal represents the population of terminals from 
presynaptic neurons onto a postsynaptic neuron. This 
induces a voltage and Ca^+response in the postsynaptic 
cell, and from the time course of ACa{t) we evaluate the 
transition rates f{t) and g{t). These enter the 'master' 
equation for the average occupations across the popula- 
tion of Ns synapses. Solving the equations for the pi{t) 
leads to our evaluation of 



which means f{t) = in our language. Here we see a 
persistent LTD dropping to 



Ns 



1 



Gampa {t) 



= goPa (t) + gipi (t) + 52P2 (t) . (22) 



LTP and LTD Induction Protocols 



1. Presynaptic Bursts 

The first protocol we used presented a burst of 
ten spikes to the presynaptic terminal and evaluated 

'^'^^jVs'^^*'' ^'^^ ^'^'^ induction period. 

The interspike interval (ISI) was constant in the burst, 
and we show in Figure 3 the value of ^"^J^ — 1 after 

the burst as a function of frequency equal to . Three 
calculations are presented. The first, shown with filled 
circles, involves the action of both the LTP inducing 
transition rate f(t) and the LTD inducing transition rate 
g{t). As in the experimental data there is a region of no 
change in AMPA conductance per synapse for very low 
frequencies, then a region of LTD until this crosses into a 
region of persistent LTP. The second calculation, shown 
with upright triangles, removes the LTD inducing tran- 
sition rate, so g{t) = 0, which is achieved by O'Connor 
et al. by the use of okadaic acid. In this calculation we 
see that LTP alone is induced at all frequencies where 
there is a measurable effect. The maximum AMPA con- 
ductance in the present model is '^'jVs'^ ~ ^ occurring 
when the lowest state is totally depleted. The value of 
unity for '^^'^'^j^^*^ — 1 is expected when a saturating 
LTP protocol is applied. Finally, a third result shown 
in Figure 3 is the set of points with inverted triangles 
which occur when one blocks kinase action, again fol- 
lowing the experimental procedures of O'Connor et al., 



-i above 

frequencies of 10 Hz. This is the smallest possible value 
in the present model, as with this induction protocol 
and f{t) — the lowest state is fully populated, and the 
AMPA conductance, in dimensionless, normalized units 
is f. 

All of this is consistent with the observations and the 
expectation of saturating LTP/LTD protocols in the dis- 
crete state plasticity model we have developed. It is im- 
portant to note that the bounded nature of the AMPA 
conductance is quite important as in many other models, 
including our own P, [TsL l25l | , there is no guarantee that 
GAMPA(t) bounded above or below. 



2. Spike Timing Plasticity 

The exploration of spike timing dependent plasticity 
at hippocampal synapses has resulted from the investi- 
gations of the phenomenon since the work of |^ 0| and 
Poo and his colleagues 0,0,01 over the past few years. 
We explored this in the present model by first present- 
ing a spike presynaptically at a time tpre and evoking 
a postsynaptic spike at tpost- The change — 1 



post 



tpre and for our model 



is a function only of r = t 
is shown in Figure 4. This reproduces the characteris- 
tic window of LTP centered near r = and of width 
« 10 ms around this point. Also shown in Figure 4 are 
the LTD regions on both sides of this window. The one 
for T negative is seen in many experiments. The LTD 
region for r positive has been seen in experiments re- 
ported by |l9j |. and it is common in models, including 
ours, which focus on postsynaptic intracellular Ca^"'"as 
inducing the chain of events leading to AMPA plasticity. 

tl9i] used cesium instead of potassium in the intra- 
cellular pipette solution, and this has been argued to 
depolarize the postsynaptic cell and broaden the action 
potential artificially l^^. To address this, has 
performed experiments in which this additional depo- 
larizing effect is mimicked by presenting a spike tim- 
ing protocol with one presynaptic spike at time tpre and 
two postsynaptic spikes with a time difference At. p9| 
uses At = 10 ms, and the outcome of this protocol for 
our model is p lotted in Figure 5 with the experimen- 
tal data ( [23; used with permission). The change in 
normalized synaptic strength resulting from this pro- 
tocol {Gampa/ Ns - 1 for the model) is shown as a 
function of the time of the second postsynaptic spike 
tpost{2) — tpre- For the experiments, normalized synap- 
tic strength is the average peak excitatory postsynap- 
tic current (EPSC) height measured 10-20 minutes after 
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the end of the pairing protocol, normahzed by the mean 
basehne peak EPSC height. It is clear that the LTP 
window is substantially larger than when we evoke just 
one postsynaptic spike and resembles the experimental 
data. We can regard this as a prediction of our discrete 
state plasticity model. Further predictions of this proto- 
col are shown in Figures 6 and 7 where At = 15 ms and 
At = 20 ms respectively. In each case there is a distinct 
LTD window for positive tpogt{2) ~ ipre and a distinctive 
dip between the LTP peaks whose separation is dictated 
by At. 

D. Synchronization of Two Periodic Neural 
Oscillators with Discrete State Synapses 

The final consequence we have investigated of our dis- 
crete state plasticity model is for the synchronization of 
oscillating neurons. We take as given that synchroniza- 
tion among populations of neurons can play an impor- 
tant role in their performing important functional activ- 
ity in biological neural networks. We have abstracted 
the synchronous activity of populations of neurons to 
the simplest setup: two periodically oscillating Hodgkin- 
Huxley (HH) neurons coupled by a synaptic current 
which we explore with and without plastic synapses. 

We have selected the postsynaptic neuron to be our 
two compartment model as described in the Appendix 
and set it into autonomous oscillations with a period 
T2 . This period is a function of the injected DC current 
into the somatic compartment. We hold this fixed while 
we inject a synaptic AMPA current 

^synapse (0 = 9AMPA{t)SA{'t){Erev — Vpost{t)), (23) 

into the postsynaptic somatic compartment. Vpost{t) is 
the membrane voltage of this postsynaptic compartment. 
gAMPA{t) is our time dependent maximal AMPA con- 
ductance, and SA{t) satisfies 

dSAjt) ^ 1 SoiVprejt)) - SAjt) ^^4) 
dt TA SiA - SoiVpreit)) 

as described in detail in the Appendix. Vprdt) is the 
periodic presynaptic voltage which we adjust by selecting 
the injected DC current into the presynaptic HH neuron. 
We call the period of this oscillation Ti . 

When gAMPA = the neurons are disconnected and 
oscillate autonomously. When gAMPAit) ^ the synap- 
tic current into the postsynaptic neuron changes its pe- 
riod of oscillation from the autonomous to T2 , which 
we evaluate for various choices of Ti. We expect from 
general arguments 10] that there will be regimes of syn- 
chronization where ^ equal integers and half-integers 



over the range of frequencies ^ presented presynapti- 
cally. This will be true both for fixed gAAiPA and when 
gAMPA varies as dictated by our model. 

In Figure 8 we present ^ as function of the frequency 
(Ti is given in milliseconds, so this is in units of 
Hz) for fixed gAMPA = 0.1^ and for gAMPA{t) = 

gAMPA^^^^f^^ determined from our model. This value 
is what we used in our earlier calculations with the two 
compartment model. It amounts to a choice for the 
baseline value of the AMPA conductance. The fixed 
gAMPA results are in filled upright triangles and, as ex- 
pected, show a regime of one-to-one synchronization over 
a range of frequencies. One also sees regions of two-to- 
one and hints of five-to-two and three-to-one synchro- 
nization. These are expected from general arguments on 
the parametric driving of a nonlinear oscillator by peri- 
odic forces. 

When we allow gAMPA to change in time according 
to the model we have discussed above, we see (unfilled 
inverted triangles) a substantial increase in the regime 
of one-to-one synchronization, the appearance of some 
instances of three-to-two synchronization, and a much 
smaller regime with two-to-one synchronization. This 
suggests that the one-to-one synchronization of oscillat- 
ing neurons, which is what one usually means by neural 
synchrony, is substantially enhanced when the synaptic 
coupling between neurons is allowed to vary by the rules 
we have described. 

We show the same results in Figure 9 for gAMPA = 
0.2^^. The fixed coupling is larger leading to stronger 
synchronization of the two neurons in a one-to-one man- 
ner even for fixed coupling. Here too (inverted, unfilled 
triangles) we see that allowing gAMPA to vary in time 
enlarges the regime of one-to-one synchronization. 

In Figures 10 and 11 we explore aspects of the internal 
dynamics of plasticity and Ca^"*" time courses for these 
results. In Figure 10 we show Ca{t) = [Ca'^~^]i{t) (scaled 
by a factor of 15 to fit on this graphic) and - 
1 in response to a presentation of periodic presynaptic 
oscillations beginning at a time 300. As noted earlier, 
the timescales for the intracellular Ca^+ processes and 
the timing in changes in GAMPA{t) are not determined 
by our model. An arbitrary constant can multiply the 
definitions of the transitions rates f{t) and g{t). Both 
quantities rapidly rise, after a small transient of LTD, to 
positive but oscillating levels. The maximum '^'^^'^ — 
1 is 1 in our model, and we see that this saturating level 
is not reached in this protocol. 

Finally, in Figure 11 we examine how the synchro- 
nization manifests itself in the postsynaptic somatic and 
dendritic compartment membrane potentials. We plot 
these potentials along with Vpre(t)- It is clear that the 
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one-to-one synchronization occurs with an in-phase os- 
cillation of the presynaptic and postsynaptic cells. The 
very short time delay between the somatic and dendritic 
compartments of the postsynaptic neuron is part of the 
model dynamics and not associated with the presenta- 
tion of periodic presynaptic spikes to the postsynaptic 
cell. The in-phase synchronization is not seen in other, 
less biophysically based, models of plastic synapses and 
represents a very desirable feature of this model. 



IV. DISCUSSION 

The observations, recent and over the years, of dis- 
crete levels for synaptic strength at individual synapses 
in the CA3-CA1 hippocampal pathways represents a fun- 
damental property important for the ways we learn and 
remember. There is a very interesting and important 
biophysical question about the mechanisms which lead to 
the expression of a few discrete levels of AMPA conduc- 
tance at which individual synapses may be found. We 
do not address this fundamental question in this paper, 
but we have used the observation to construct a model 
based on discrete levels with transition rates among the 
levels determined by biophysical dynamics. 

We have formulated a discrete level synaptic system 
in a general way with L levels allowed to the AMPA 
conductance, and then, following the observations of 
O'Connor et al., we specialized to L = 3. The dynamical 
variables in our model when L — Z are the average occu- 
pation numbers of each level 'Pit) — [po(0:Pi(^)i-P2(i)]- 
These are averages over a collection of Ns synapses 
which contribute to the overall AMPA determined re- 
sponse of the neuron. While each individual synapse 
resides in one of three discrete states, so the individual 
occupation numbers at any given synapse are either zero 
or one, the average occupation numbers are smoothly 
varying, subject only to pn{t) + pi{t) + P2{t) = 1, by 
definition. 

We developed differential equations for P{t) which are 
linear in the pi{t), 1 — 0,1,2 



dpijt) 
dt 



Y^Mu-pvit), 



(25) 



;'=o 



and where the transition rates Mui are determined 
by nonlinear membrane voltage and intracellular 
Ca^+ dynamics. 

From the observations of O'Connor et al. we argued 
that the transition rates shown in Figure 1 sufficed to 
explain their measurements, and using their reported re- 
sults we were able to determine that the conductances of 
the three individual levels in normalized, dimensionless 



units were go = f 1 5i = 2, 52 = 2. The time dependence 
of the normalized, dimensionless AMPA conductance per 
synapse is then 



Gam PA {t) 



(26) 



Using our values for the gi and a dynamical model of 
the transition rates f{t),g{t) as shown in Figure 1, we 
reproduced the observed plasticity in response to a burst 
of presynaptic spikes with interspike intervals (ISIs) over 
the observed range. Further we made predictions for the 
response of this model to spike timing plasticity both 
for one presynaptic and one postsynaptic spike and for 
the case of two postsynaptic spikes evoked At apart ac- 
companied by one presynaptic spike. We presented our 
resuhs for At = 10, 15, ...,20 ms. 

Finally we examined the dynamical role played by 
this discrete state plasticity model in the synchronization 
of two periodically oscillating Hodgkin- Huxley neurons. 
One such neuron oscillating with period was driven by 
another such neuron with period Ti . The final period T2 
of the driven neuron, relative to Ti , was plotted against 
^ and showed familiar regions of synchronization. For 



fixed AMPA coupling qampa 



Ns 



we found syn- 



chronization over some range of ^ and then demon- 
strated that allowing qampa to vary according to the 
plasticity model resulted in a much larger regime of one- 
to-one synchronization with the two neurons oscillating 
in-phase. The results for synchronization have not been 
tested experimentally, though some experiments using 
dynamic clamp based synapses have been performed. 

One striking aspect of the discrete state model, cer- 
tainly not limited to our own work, is that the AMPA 
conductance has natural upper and lower bounds. Many 
other models of plasticity, including our own, do not 
share this important feature. 

Some of our results, in particular the strengths of the 
normalized, dimensionless conductances of the synaptic 
levels are dependent primarily on the data of O'Connor 
et al. All of the transition rates are determined by our 
two compartment model for the neuron, as presented in 
the text and in the Appendix. 

Our model exhibits a number of features of LTP and 
LTD observed experimentally at CA3-CA1 synapses, in- 
cluding trapping of synapses in a high-strength state, 
separability of potentiation and depression by simulated 
inhibition of kinase or phosphatase activity, and spike 
timing-dependent plasticity (O'Connor et al.; |29|). 
However, there are a number of ways in which the model 
can be developed further. For example, using the pro- 
tocols of O'Connor et al., population LTP at CA3-CA1 
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synapses rises gradually to a peak level over a few min- 
utes; LTD takes a few minutes longer than this to de- 
velop fully. Our model does not yet include such a long 
timescale, but could be modified phenomenologically to 
do so. O'Connor et al. have also found that the high 
locked-in state in populations of synapses builds up over 
several minutes. To model this phenomenon, we would 
again need to include a longer timescale. 

While the spike timing-dependent plasticity induced 
in our model by 1 presynaptic and 2 postsynaptic spikes 
is similar to that observed by |29j , our result for 1 presy- 
naptic and 1 postsynaptic spike appears to differ from ex- 
perimental observations (G. M. Wittenberg and S. S.-H. 
Wang, unpublished data used with the authors' permis- 
sion). In particular, they observed little LTP but signif- 
icant LTD near tpost - tpre = 0. At other values of tpost - 
tpre, they observed only LTD. This suggests that such a 
protocol provides insufficient postsynaptic Ca^+influx to 
induce LTP reliably. In contrast, our model shows a nar- 
row but clear window of LTP centered near tpost - tpre = 
0. At these values of tpost - tpre, the Ca^+influx in our 
model is sufficient to give f a relatively large value and 
thus induce LTP. If the data of Wittenberg and Wang 
are correct, then our model will need to be adjusted so 
that this Ca^+influx is not sufficient to induce LTP. 

In addition, the model can be refined to match the 
results of various LTP and LTD induction protocols 
that we have not simulated here but that are used by 
O'Connor et al. and others in their experiments, like 
theta burst stimulation and pairing protocols. In the 
long term, a model that more accurately describes the 
postsynaptic signaling pathways will eventually account 
for all of these various features of the data in a biologi- 
cally satisfying manner. 

While our model is based on the idea of discrete state 
synapses, by design it describes only populations of such 
synapses. In future work we plan to address this discrete- 
ness at the level of single synapses and small numbers of 
synapses. In a single-synapse model, transition proba- 
bility functions would likely replace the transition rate 
functions we have described in the present model, but 
many other features of the model could be retained. 

This model will change over time and be improved by 
further understanding of the biophysical processes lead- 
ing to the discrete states and their transitions among 
themselves. The general framework we have presented 
describing how the three observed states are connected 
and several general results about that system will remain 
as the representation of the transition rates is improved. 
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V. APPENDIX 

In this Appendix we give the details of the two 
compartment model of the postsynaptic neuron used 

in our calculations. The somatic compartment is the 
site of spike generation by the familiar Hodgkin-Huxley 
(HH) Na+, K+, and leak currents. The dendritic 
spine compartment has these currents as well as a 
voltage gated calcium current and glutamate driven 
NMDA and AMPA channels. The Ca^+dynamics in 
the dendritic compartment drives the synaptic plastic- 
ity, namely changes in the AMPA conductance. 

A. The Somatic Compartment 

The dynamical equation for the somatic compartment 
takes the general form, 

Cm^^ = lNa{Vs{t),t)+lKiVs{t),t)+lL{t) 
+ Isdc + Is{t) 

+ Gs^D{VDit)-Vsit)) (27) 

The currents Inu, Ik, II, are the familiar HH Na"*", K+, 
and leak currents. Isdc is a DC current used to set the 
resting potential of the cell. Is{t) is an externally man- 
aged time dependent current injected into the somatic 
compartment. It allows us to evoke an action poten- 
tial at a specific time in the somatic compartment. This 
propagates back to the dendritic compartment to induce 
a depolarizing effect. Gs^DiVoit) — Vs{t)) represents 
the current flowing into the somatic compartment from 
the dendritic compartment. It couples the voltages of the 
somatic and dendritic compartments. Cm is the mem- 
brane capacitance. 

The value of the currents I^a, Ik, and II are deter- 
mined as usual with 

lL{t) = 9l{El - Vs{t)), (28) 

where is the conductance of the leak current and El 
is the reversal potential. The voltage gated currents are 
described by 

I{V,t) = -gg{t,V){E,^-V), (29) 

where E(,q is the reversal potential and g is the maximal 
conductance. Both these values are fixed. The value of 
g{t, V), the fraction of open channels, on the other hand 
depends on the membrane potential and time. 

In the case of channels in which g(t, V) changes, the 
value of g(t, V) depends on the state of 'gating parti- 
cles' m{t, V) and h{t, V), where m{t, V) is the activation 



gate and h{t, V) represents the inactivation gate. If A'' 
is the number of activation gates, and M, the number of 
inactivation gates then 

g{t,V) = m{t,V)''h{t,V)'' 

The state of the gating particles, given by m{t, V) or 
h{t, V), is a function of the membrane potential as well 
as time. 

These gating variables, denoted by X{t), are taken to 
satisfy first order kinetics. 

dX{t) ^ Xo{V{t))-X{t) 
dt Tx{V{t)) 

= axiVmi - X{t)) - Px{V{t))X{t){30) 

From the standard HH model, we have the following rela- 
tions for the conductances of the Na'^ and currents. 

gNa{V,t) = m{V,tfh{V,t) 
gK{V,t) = n{V,tf 

where ■m{V, t),n{V, t) are activation gating particles and 
h{y,t) represents the inactivation gating particle. At 
the end of this appendix we give the functions Xq.tx, 
or ax{V)^ Px{y) for each of the voltage gated ionic cur- 
rents, in addition to listing all the model parameters used 
in the simulations. 



B. The Dendritic CompEirtment 

The dynamics of the dendritic compartment mem- 
brane potential is given by 

Cm^^^ = lNa{VD{t),t) + lK{VD{t),t) + iLiVoit)) 

+ Ia{Vd it) , t) + Im ( Vb [t] , t) + iDdc 
+ lAMPA{t) -I- lNMDA{t) + IvGCc{t) 

+ GD^s{Vs{t) - VdH)) (31) 

lodc is a DC current used to set the resting potential 
of the cell. In a, Ik and 1^ represent the standard HH 
ionic and leak currents used in the somatic compartment 
as described above. In addition to these ionic currents 
we have considered two additional currents, Ia and 
Im- I a currents have been reported to modulate the 
width of action potentials and influence the excitability 
of the cell. In our model, I a attenuates the dendritic 
action potential, which is evoked by backpropagation of 
the somatic action potential. 

The gating equations for Im and I a are, 

gM{t, V)=u{t,Vf 
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gA{t,V) ^ a{t,V)b{t,V) 

where u(t, V) and a(t, V) are activation gating particles 
and b{t, V) represents an inactivation gating particle. 

As mentioned above all ionic currents in the den- 
dritic compartment are also given in terms of Ohm's 
law, (123 ■ Again G£i^s(Vs(i) — ^oit)) represents cur- 
rent flowing into the dendritic compartment through the 
somatic compartment. 

In addition to the currents mentioned above we have 
three other currents critical to the synaptic plasticity 
discussed here. There is a current associated with the 
ligand gated NMDA receptors (NMDARs). The form 
for this is 

iNMDAit) = gNMDASN{t)B{VD{t)){VNMDA-eq " VoitW 

where gNMDA is the maximal conductance associated 
with the channel. S]\[{t) ranges between zero and unity, 
representing the percentage of open channels at any 
time. To achieve the time course of this process in NM- 
DARs, we use a two component form for Sn, 



This voltage dependent conductance depends on the 
extracellular magnesium concentration. The voltage de- 
pendence of the current is mediated by the magnesium 
ion which, under normal conditions, blocks the channel. 
The cell must therefore be sufficiently depolarized to re- 
move the magnesium block. Finally for this excitatory 
channel Vjvmda-gij ~ OmV. 

I AM PA represents the ligand gated AMPA receptor 
current. This is taken to be of the form, 

IaMPA — gAMPASA{t)(yAMPA~eq - ^d(O) (36) 

where gAMPA is the maximal conductance for this chan- 
nel and SAit) is the percentage of open channels, satis- 
fying 



S'Ar(i) = WfSNlit) + (1 - Wf)SN2{t), 

< Wf < 1, and where SNi{t), 1 — 1,2 satisfies 
dSm{t) _ 1 SaiVpreit)) - Smit) 



dt 



TNI SiNl — So{Vpre{t)) 



(33) 



(34) 



Vpre is scaled to lie between and 1 as it represents the 
arrival of an action potential at the presynaptic terminal. 
Its function is in turn to release neurotransmitter. 

SoiVpreit)) is a step function which rises sharply from 
to 1 when neurotransmitter is released as a result 
of the presynaptic action potential. When this occurs 
S'jv;(i) rises from zero towards unity with a time con- 
stant TNi{SiNi — 1). When the effect of presynaptic ac- 
tion is completed, SNi(t) relaxes towards zero with a 
time constant tniSini- Wf represents the fraction of 
fast NMDA component contribution to NMDA current. 
In our model we have chosen Wf — .81, tni — 67.5, 
SiNi = 70/67.5, TN2 = 245, Sin2 = 250/245. In addi- 
tion the conductance of the NMDA current depends on 
postsynaptic voltage via the term B{V) whose form is 
given as. 



dSAjt) ^ 1 S„{Vpre{t)) ~ SAjt) 
dt TA SiA - SoiVpreit)) 



(37) 



Again the rise time is less than a millisecond. In our 
formulation this time is ta{Sia ~ 1), which we set to 0.1 
ms. AMPA currents decay in approximately 1-3 ms. In 
our formulation this decay time is taSia, which we set 
to 1.5 ms. We also take VAMPA-eq — mV. 

The final and very important ingredient in inducing 
synaptic plasticity is the voltage gated calcium channel 
(VGCC) . We have used the low threshold current It for 
this. The current from this channel takes the form, 

IvGCcit) = gcGHK(V{t))ml{t)h,{t) (38) 

where gc is the maximal conductance of this channel, 
mc(t) is the activation function, and hc{t) is the inacti- 
vation function. GHK(V) is the Goldman-Hodgkin-Katz 
function, 



GHK{V) 



V [Ca^+] 



it) 



e RT 



Co 1- 
V Ca{t) - [Ca^ 



-2VF 



Cn 



1 



-2VF 



(39) 



which is used because of the disparity in the intracellu- 
lar [Ca^+] i and the extracellular [Ca^+]o concentrations. 
F is Faraday's constant, R is the gas constant, and T 
the absolute temperature. Other factors of the GHK 
equation are absorbed in the conductance gc- Co is the 
equilibrium intracellular [Ca^+] concentration, which is 
about 100 nM. 



B{V) = 



.288[Mg2+]e- 



.062V ' 



(35) 



where the concentration of magnesium is in mM and the 
voltage is in mV. For simulation purposes we have taken 
the physiologically reasonable value of [A/g^+] = 1 mM. 



C. Coupling between the somatic and dendritic 
compartments 

The coupling parameters between the two compart- 
ments are determined from the cytoplasmic resistance 
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and the metric dimensions of each compartment. We 
take the specific cytoplasmic resistance of the cell to be 
n = 200f^ cm. 

We take the somatic compartment to be a isopotential 
sphere of dsoma = 32.5^m in diameter and the dendritic 
compartment to be an isopotential cylinder of diameter 

doendrite = lOyUm and length loendrite = 360/XTO. 

In order to determine the coupling resistance value 
we assume the somatic compartment to be a cylinder of 

equivalent surface area. We then have the total cyto- 
plasmic resistance of the somatic compartment, 



C a AM PA = gAcSA{t){VNMDA- 



Riso 



Tilso 



Crosssection( Asoma ) 



.007839 X 



while the total cytoplasmic resistance for the dendritic 
compartment is 



R 



Til Dendrite 



IDendrite 



1.713 X lo'^n 



Crosssection( A/Jen ) 
Therefore, the average cytoplasmic coupling resistance 

n Rl Soma ~^ Rl Dendrite Rl Dendrite o/^i tr\7r\ 

Rj = « = .861 X 10 s2 

2 2 

The coupling parameters in units of mS/cm? are calcu- 
lated as 



and 



Gs^D = 



Gd^s 



1 



■^SomaRl 



= 3.5mS/cTn^ 



^Den 



Ri 



= l.OmS/cm^ 



(40) 



(41) 



D. Calcium Dynamics 

The dynamics of intracellular calcium [Ca^^(i)] in 
the dendritic compartment, which affects the efficacy of 
synaptic strength, is comprised of [Ca'^'^{t)] decaying to 
an equilibrium value of Co on a timescale of tc ~ 15ms, 
which we take to be about 30 ms in our model, plus fluxes 
of [Ca^+{t)] due to the three channels, AMPA, NMDA, 
and VGCC considered in the dendrite model above. The 
first order differential equation for [Ca'^~^{t)] = Ca{t) 
then is 



dCajt) 
dt 



— {Co 

TC 

GavGCC 



Ca{t)) + Cgnmda + Cgampa 

(42) 



where 



eg 



Vdenit)) 



CavGCC = gccGHK{V{t))ml{t)K{t) 

The constants QnCi 9ac, gcc are not the same, even di- 
mensionally, as the conductances in the voltage equa- 
tion. Their values, given at the end of this appendix, 
reflect among other things, that the net AMPA current 
is composed primarily of other ions besides Ca^+ and 
that NMDA channels are highly permeable to Ca^+ ions. 
This completes the description of our model. 



E. Parameters in the Two Compartment Model 

The various constants appearing in our two compart- 
ment model are collected here. 

• The somatic compartment : The maximal con- 
ductances, in units of mS/crn^, of the ionic cur- 
rents are, qnu = 215, = 43, — .813. The 

reversal potentials in units of mV are Vjsia-eq = 



50, Vk- 



^ -95, and Vl-co = -64. The DC 



CttNMDA = 9NcSNit)B{yden{t)){yNMDA-eq — Vdenit)) 



current injected into the somatic compartment 
Isdc = — 7.0/Lt^/cm^, so that the cell is at -75mV 
at rest. The magnitude of the additional current 
injected into the somatic compartment is Isoma = 
160.8/xA/cm^; in using it to induce a postsynaptic 
spike, it is taken to have a duration of 1 ms or less. 

• The dendritic compartment : For the standard 

HH ionic currents we have the same parameters 
as above. Maximal conductances associated with 
various dendritic currents, in units of mS/crn?, are 
9NMDA = 0.05, gAMPA = 1-75 and ,gc = 1- x 10"^. 
In the Mg^^ blockage fimction B{V), wc take 
= imM. In the GHK function, the ra- 
tio of external [Ca^+] to equilibrium intracellular 
[Ca?+] is 15000. The temperature is 25°C. The 
conductance values for additional potassium cur- 
rents used are gu = 6.7 and gA = 100 in units of 
mSjcrr?. Finally, ludc = -7.0/xA/cm^ 

• Calcium dynamics : For calcium dynamics we 

have Tc = 30ms, gj\[c — . 15, . g^c = 1-5 x 10^"'', 
and gcc = 3.5 x 10~^. These are in units of 
mV~^ms~^. Co, the basal calcium concentration 
in the cell, is normalized to 1. 

Finally in the equation for the NMDA and AMPA chan- 
nel open percentages S'jv(i) and SA{t), we use the 'step' 
function 

So{V) = i(l + tanh(120(y - .1))) 
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and the constants ta = lAms, TNfast = 67.5ms, 

TNsIow = 245.0ms, Sia = ' ^^Nfast = 

c _ 250 

^INsloW — OAK • 



F. Activation and deactivation parameters of 
various channels 

f^-. _ .32{13-{V-Vth)) _ 0.28((V-yth)-40) 

C-myV J — (iz-{v--vtWT) Pm\V ) — ((v-vth)-4.o) 

e 4.0 -1 e 5 -1 

ah{V) = .1286 """^s"*"' ^h{V) = 4o-(v-v.,o 



an{V ) — — (i5-(v-vthT) Pn\V ) 



(15-(V-Vtfe)) Pnyv ) — - (V-Vtfe)-lO 

e 5 — 1 e 40 

'TT'Co(l^) = ("52+ V) Tmciy) = .204 H (i5+v) 

1+e 6.2 e 16.7 +e is. 2 

/ico(l^) = (72 + V) Q!u(l^) = (V+5^2.7r 

1+e 4 e 23 

« rv^ - ,0 016 ^ _ -0-05(V;+20) 

e — — e 15 —1 

Pa\y ) — (v+10) OibyV ) — (v + is) 

e 8 — 1 e 15 

PbiV) = (yl^js) 

e" 12 +1 

(V + 466) 

Thc{V) = 0.333e-ss:^ if F < -81 

= 9.32 + 0.333e-^T^ iiV > -81 



where Vth = -48 mV. 
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VI. FIGURES 




FIG. 1: Three state model of synaptic plasticity. Individ- 
ual synapses can move among the three states, marked for 
the "low" state, 1 for the "high" state, and 2 for the "high 
locked-in" state. The rules for state transition depend on the 
transition rates, f{t) and g{t), governed by changes in in- 
tracellular calcium concentration. These transition rates are 
determined by LTP and LTD induction protocols. 
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FIG. 2: Steady state values for the transition rates / and 
g plotted as function of the magnitude of the change in in- 
tracellular calcium concentration. ACa is in arbitrary units. 
For small changes in intracellular calcium concentration only 
g is nonzero, corresponding to LTD induction, and for larger 
changes in intracellular calcium concentration, both / and g 
are nonzero with / being greater than g, corresponding to an 
LTP induction protocol. 
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FIG. 3: Frequency-plasticity curve. The change in normal- 
ized AMPA conductance per synapse, Gampa/Ns — 1, is 
plotted as function of the frequency of a periodic burst of 
10 presynaptic spikes presented to the presynaptic termi- 
nal. The circles represent synaptic plasticity for the full three 
state model. The upward-pointing triangles represent synap- 
tic plasticity with the term g{t) set to 0, corresponding to 
blocking phosphatase activity in the postsynaptic cell. One 
sees and expects LTP alone. The downward-pointing trian- 
gles represent the change in synaptic plasticity with the term 
f{t) set to 0, corresponding to blocking kinase activity in the 
postsynaptic cell. We observe and expect LTD alone in this 
case. These results are quite similar to the observations of 
O'Connor et al. 




FIG. 4: Spike timing dependent plasticity protocol. The 
change in normalized AMPA conductance per synapse, 
Gampa/Ns — 1, plotted as a function of the delay, r = 
tpoat — tpre (ms) , between presentation of a single presynaptic 
spike at tpre and postsynaptic spike at tpost- 
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FIG. 5: Change in normalized synaptic strength as a func- 
tion of the delay t = tpost(2) — tpre when a single presynaptic 
spike is paired with two postsynaptic spikes 10 ms apart. 
tpost(2) is the time of the second postsynaptic spike. Model 
results {points connected by lines) are plotted with exper- 
imental data of G. M. Wittenberg (large filled circles with 
error bars; Wittenberg 2003, used with permission). Nor- 
malized synaptic strength, for the model, is the normalized 
AMPA conductance per synapse (Gampa/Ns) after the pair- 
ing. For the experiments, it is the average peak excitatory 
postsynaptic current (EPSC) height measured 10-20 minutes 
after the end of the pairing protocol, normalized by the mean 
baseline peak EPSC height {Error bars: standard error of the 
mean). In the experiments, pairing was repeated 100 times 
at 5 Hz. 
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FIG. 6: The change in normalized AMPA conductance per 
synapse, Gampa/Ns — 1, plotted as a function of the delay 
T = tpost(2) ^ tpre (uis) wlicu wc prcscnt a single presynaptic 
spike and a pair of postsynaptic spikes with a fixed time dif- 
ference of 15 ms. tpost(2) represents the time of presentation 
of the second postsynaptic spike. A distinct dip in the po- 
tentiated AMPA conductance is observed for times when the 
presynaptic spike falls in the time between the 2 postsynaptic 
spikes. 
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FIG. 7: The change in normaUzed AMPA conductance per 

synapse, Gampa/Ns ~ 1, plotted as a function of the delay 
T = tpost(2) ~ tpre (his) wheu WO sclcct a single presynaptic 
spike and a pair of postsynaptic spikes with a fixed time dif- 
ference of 20 ms. tpost{2) represents the time of presentation 
of the second postsynaptic spike. A distinct dip in the po- 
tentiated AMPA conductance is observed for times when the 
presynaptic spike falls in the time between the 2 postsynaptic 
spikes. 
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FIG. 8: T1/T2, the ratio of the interspike interval Ti of the 
presynaptic neuron to the interspike interval T2 of the post- 
synaptic neuron, is plotted as a function of the presynaptic 
input frequency, 1000/Ti Hz, for a synapse starting at a base 
AMPA conductance of g amp Ait = 0) =0.1 mS/cm^. We 
see that the one-to-one synchronization window is broadened 
when the static synapse is replaced by a plastic synapse as 
determined by the three state model. 
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FIG. 9: T1/T2, the ratio of the interspike interval Ti of the 
presynaptic neuron to the interspike interval T2 ms of the 
postsynaptic neuron, is plotted as a function of the presy- 
naptic input frequency, 1000/Ti Hz, for a synapse starting at 
a base AMPA conductance of gAMPA{t = 0) = 0.2 mS/cm? . 
We see that the one-to-one synchronization window is broad- 
ened when the static synapse is replaced by a plastic synapse 
as determined by the three state model. 



FIG. 11: Vsoma(t), ^dendrite (*) ^^'^ ^^pre(^), plotted as func- 
tions of time, when the presynaptic and postsynaptic neu- 
rons are synchronized. Note that the presynaptic and post- 
synaptic neurons are synchronized in-phase with an internal, 
Vsoma(t) to V'gpjjjg(f), time difference determined by the two 
compartments of the model neuron. 
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FIG. 10: Intracellular calcium concentration, scaled by 
15, and the change in normalized synaptic strength, 
GAMPA{t)/Ns — 1, is plotted as a function of time in the 
case when a periodically spiking postsynaptic cell is driven 
by a periodically spiking presynaptic input. 
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